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Abstract 

Streambed scour at bridge piers is among the leading 
causes of bridge failure in the United States. Several pier-scour 
equations have been developed to calculate potential scour 
depths at existing and proposed bridges. Because many pier- 
scour equations are based on data from laboratory flumes and 
from cohesionless silt- and sand-bottomed streams, they tend to 
overestimate scour for piers in coarse-bed materials. Several 
equations have been developed to incorporate the mitigating 
effects of large particle sizes on pier scour, but further investi- 
gations are needed to evaluate how accurately pier-scour depths 
calculated by these equations match measured field data. 

This report, prepared in cooperation with the Montana 
Department of Transportation, describes the evaluation of five 
pier-scour equations for coarse-bed streams. Pier-scour and 
associated bridge-geometry, bed-material, and streamflow- 
measurement data at bridges over coarse-bed streams in Mon- 
tana, Alaska, Maryland, Ohio, and Virginia were selected from 
the Bridge Scour Data Management System. Pier scour calcu- 
lated using the Simplified Chinese equation, the Froehlich 
equation, the Froehlich design equation, the HEC-18/Jones 
equation and the HEC-18/Mueller equation for flood events 
with approximate recurrence intervals of less than 2 to 100 
years were compared to 42 pier-scour measurements. Compar- 
ison of results showed that pier-scour depths calculated with the 
HEC-18/Mueller equation were seldom smaller than measured 
pier-scour depths. In addition, pier-scour depths calculated 
using the HEC-18/Mueller equation were closer to measured 
scour than for the other equations that did not underestimate 
pier scour. However, more data are needed from coarse-bed 
streams and from less frequent flood events to further evaluate 
pier-scour equations. 



Introduction 

Streambed scour at bridge piers is among the leading 
causes of bridge failure in the United States (Landers and Muel- 
ler, 1996, p.l). As a result, the Federal Highway Administra- 
tion (FHWA) has developed methods for State highway agen- 



cies to calculate potential scour depths at existing and proposed 
bridges (Richardson and Davis, 2001). Methods for calculating 
pier scour are based on empirical equations relating maximum 
scour depth to various hydraulic and bridge-geometry variables. 
Many of the empirical equations currently available (2004) are 
based on data from laboratory flumes and from cohesionless 
silt- and sand-bottom streams. Equation improvements con- 
tinue to be made with the ultimate goal of minimizing under- 
estimation and overestimation of scour. If pier scour is under- 
estimated, scour depths assumed in the bridge design could be 
exceeded during large floods. Excess scour could lessen sup- 
port for the bridge pier and destabilize the bridge. To ensure 
that pier-scour depth is not underestimated, some empirical 
equations have been adjusted to yield more conservative 
(larger) scour estimates. However, calculated scour depths 
from those equations might indicate that more bridges are 
scour-critical (subject to failure due to scour) than is actually the 
case and thus may lead to expensive over-design or unnecessary 
retrofitting of pier foundations. The goal for design of new 
bridges and analysis of existing structures is to ensure that 
bridge foundations withstand the effects of scour, but are not 
larger, deeper, or more expensive than necessary. Therefore, 
estimated scour depths from a pier-scour equation for a given 
set of site and flood conditions need to be as accurate as possi- 
ble. However, when estimates are in error, scour needs to be 
overestimated rather than underestimated for safety consider- 
ations in bridge design. 

To improve the understanding of scour processes and to 
develop more reliable pier-scour equations, the FHWA and 
many State highway agencies have cooperated with the U.S. 
Geological Survey (USGS) to collect onsite scour data at 
bridges. Data from these studies were analyzed by Landers and 
Mueller (1996) and used to evaluate 14 empirical pier-scour 
equations. Three hundred eighty -four pier-scour measurements 
at 56 bridges in 13 states were then compared to pier scour cal- 
culated by each of the 14 equations. Comparisons showed that 
none of the 14 equations accurately calculated scour for the full 
range of conditions measured in the field. Moreover, Landers 
and Mueller (1996, p. 111-112) found that the pier-scour equa- 
tion recommended in the FHWA Hydrologic Engineering Cir- 
cular 18, Second Edition (Richardson and others, 1993) overes- 
timated pier scour for many measurements. This equation. 
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referred to as the HEC-18 (2nd edition) equation in this report, 
did not account for the effects that coarse streambed material 
has on scour. 

For coarse-bed streams. Landers and Mueller (1996, p. 95- 
119) identified several empirical equations that accounted for 
bed-material size, including the Simplified Chinese equation, 
the Froehlich equation, and the Froehlich design equation. 
Landers and Mueller (1996) cite Gao and others (1993) as the 
source for the Simplified Chinese equation and Froehlich 
(1988) as the source for the Froehlich and Froehlich design 
equations. In the Landers and Mueller study (1996, p. 109- 
111), the Simplified Chinese equation and the Froehlich equa- 
tion frequently underestimated pier scour, while the Froehlich 
design equation only rarely underestimated pier scour. 

Since publication of the HEC-18 (2nd edition) (Richard- 
son and others, 1993), a third edition (Richardson and Davis, 
1995) and a fourth edition (Richardson and Davis, 2001) have 
been published. The HEC-18 equation in the third edition 
included a pier-scour correction factor (^"4) for coarse bed mate- 
rial developed by J. Sterling Jones, Federal Highway Adminis- 
tration. The fourth edition included a modified K^ factor that 
better accounted for streambed armoring in coarse-bed chan- 
nels. Application of this modified ^4 factor was expected to 
result in substantially smaller and more reliable estimates of 
pier-scour depths in coarse-bed streams (Mueller, 1996). 

The third- and fourth -edition versions of the HEC-18 
equation are hereinafter referred to in this report as the HEC-18/ 
Jones equation and the HEC- 1 8/Mueller equation, respectively. 
Because of a small difference in the equation for critical veloc- 
ity, the HEC-18/Mueller equation used in this report (Mueller, 
1996, p. 158-160) is slightly different from the HEC-18/Muel- 
ler equation included in Richardson and Davis (2001, p. 6.6). 

Because they account for bed-material size, and in 
some instances streambed armoring, the Simplified 
Chinese, Froehlich, Froehlich design, HEC-18/Jones, and 
HEC-18/Mueller equations are generally considered appropri- 
ate for calculating pier scour in coarse-bed streams in mountain 
states like Montana. However, the five equations have not been 
evaluated for use in coarse-bed streams. Therefore, the USGS 
and the Montana Department of Transportation (MDT) initiated 
a cooperative study in July 2000 to evaluate the five pier-scour 
equations by comparing calculated pier scour with measured 
pier scour from coarse-bed streams. 



were evaluated based on the accuracy of pier-scour estimates 
and number and magnitude of underestimates resulting from 
each equation. 

The definition of coarse-bed streams has evolved. Rich- 
ardson and Davis (1995, p. 38) indicated that the pier-scour cor- 
rection factor for coarse-bed streams (^4) should be used only 
where D^q (the particle size for which 50 percent is finer) is 
greater than 60 mm. Later, Richardson and Davis (2001) char- 
acterized coarse-bed streams as those where D^q is greater than 
40 mm. They further indicated that application of the K^ factor 
was appropriate where D50 is greater than 2 mm and ZJgj is 
greater than 20 mm. In this report, coarse-bed streams were 
considered to be those having Djg greater than 50 mm. Thus, 
all sites analyzed easily met the more recent criterion for coarse- 
bed streams given by Richardson and Davis (2001). 



Description of Pier-Scour Data Used to 
Evaluate Equations 

Coarse-bed (D^q > 50 mm) pier-scour data compiled for 
this study were selected from the USGS Bridge Scour Data 
Management System (BSDMS) (Chad R. Wagner, U.S. Geo- 
logical Survey, written commun., 2002). The selected BSDMS 
data included 15 pier-scour measurements from two Montana 
sites: Gallatin River at U.S. 191 near Gallatin Gateway (site 
32), and Yellowstone River at U.S. 89, near Emigrant (site 33) 
(tables 1 and 2). The data selected from the BSDMS also 
included measurements (not previously analyzed by other 
researchers) made at site 33 during 1996 and 1997, when both 
peak discharges were close to the 50-year and 100-year floods, 
respectively. Furthermore, scour data for the Bitterroot River 
at U.S. 93 near Darby, Mont, (site 85), were recently added to 
the BSDMS and are included in this study. 

This study also used selected BSDMS data from four other 
States. Included were 10 measurements from 2 rivers in 
Alaska — the Susitna River (site 1) and the Tazlina River (site 
4); 4 measurements from the Youghiogheny River in Maryland 
(site 23); 1 measurement from the Little Miami River in Ohio 
(site 44); and 9 measurements from 2 rivers in Virginia — the 
Tye River (site 53) and Reed Creek (site 55). 



Purpose and Scope 



Methods of Data Collection 



This report describes the results of an evaluation of five 
pier-scour equations for coarse-bed streams. Scour depths cal- 
culated by the equations were compared with 18 pier-scour 
measurements at 3 bridge sites in Montana, 10 measurements at 
2 sites in Alaska, 4 measurements at 1 site in Maryland, 1 mea- 
surement in Ohio, and 9 measurements at 2 sites in Virginia 
(fig. 1 and table 1). Scour measurements for each site are sum- 
marized in tables 1 and 2. The Simplified Chinese equation, the 
Froehlich equation, the Froehlich design equation, the 
HEC-18/Jones equation, and the HEC-18/Mueller equation 



Scour measurements in Montana were made using stan- 
dard USGS sounding equipment that included either a four- 
wheel base or a bridgeboard device with a cable-suspended 
Columbus sounding weight attached to an A- or B-type sound- 
ing reel deployed from the upstream side of each bridge. Initial 
soundings typically were made for the entire cross section at the 
upstream face of the bridge before the runoff period to obtain 
baseline data. Velocity measurements were made using a ver- 
tical-axis current meter (Rantz, 1982) at several stations on both 
sides of each pier, outside the zone where flow typically 
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STATE WITH PIER-SCOUR DATA USED 
IN THIS STUDY AND NUMBER OF 
SITES ANALYZED 



Figure 1. States from which pier-scour data were compiled. 



accelerates near the pier. Sounding depths were related to a sur- 
veyed water-surface elevation, which was then related to a 
nearby vertical datum. Scour and velocity measurements at 
sites in the other States were made in a similar manner. Once 
the baseline cross-section data were plotted, the reference sur- 
face (Landers and Mueller, 1993, p. 2075-2080) was estimated 
and used with the scour soundings to derive the estimated scour 
depth. When more than one high-discharge scour measurement 
was made at a bridge, each was considered to be an independent 
measurement. 

Bed material for each Montana site was characterized by 
particle counts (Wolman, 1954) performed near the bridge 
piers. Though the precision involved in determining bed-mate- 
rial size depends on a number of factors related to size distribu- 
tion, spatial variation, sample size, and technique (Wolman, 
1954; Kellerhals and Bray, 1971; Hey and Thorne, 1983), val- 
ues are typically reported to two or three significant figures. 
Methods used to estimate bed-material sizes at sites outside of 



Montana varied, but are believed to be consistent with those 
applied in Montana. 

Comparison of Data Sets Used to Develop and 
Evaluate Equations 

Though researchers in several States are working to col- 
lect, compile, and analyze bridge-scour data, the database con- 
taining pier-scour measurements and associated site and flow 
information for coarse-bed streams is fairly small. Therefore, 
an evaluation database that is totally independent of the data- 
base used to develop the pier-scour equations was unavailable 
for this study. In this study, 17 of the 42 measurements used to 
evaluate pier-scour equations also were used to develop the 
HEC-18/Mueller equation (D.S. Mueller, U.S. Geological Sur- 
vey, written commun., 2002). However, 25 of the 42 pier-scour 
measurements used in this evaluation were not used in the 



Table 1. Hydrologic and hydraulic data for sslsctad pier-scour sitos 



[Drainage area and channel-slope data for site 85 from Bridge Scour Database Management System. Area and slope data for other bridges from Landers and Mueller (1996). Discharge measurement 

data from Bridge Scour Database Management System, unless otherwise noted. Abbreviations: BSDMS, Bridge Scour Database Management System; mi , square miles; ft/ft, foot per foot; 

ft /s, cubic foot per second; NA, not applicable; Q.,, approximate peak discharge at the x-year recurrence interval, in ft /s; S.R., State Road; U.S., United States; USGS, U.S. Geological Survey. Symbol: <, less 

than; ~, data not available] 
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Table 2. Summary of pier-scour measurements 

[Pier identification: number or location of pier at which measurement took place. Scour condition: refers to general condition of the bed-material movement upstream from the bridge 
during the time of measurement — clear, bed material not in motion; live, bed material generally in motion; *, scour condition based on comparison of measured velocity and computed 
incipient-motion velocity, rather than observation or hydrologist's judgement. Pier shape: sharp-nosed, round-nosed. Pier skew: measured in degrees from line parallel to flow. Measurement 
error: scour measurement error (in feet) estimated by hydrologist. For example, for measurement number 228 the hydrologist measured 0.8 ft of scour and estimated a measurement error of 
0.3 ft. Therefore, the measured scour depth could vary from 0.5 ft to 1.1 ft. Abbreviations: no., number; L, length; ft, foot; deg, degrees; b, width; V, velocity; ft/s, foot per second; mm, 
millimeters. Symbol: — , not available] 
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Clear 


Sharp 


20.0 





5.0 


11.5 


17.0 


70 


90 


96 


2.0 


1.0 




7 


07/02/71 


P4 


Clear 


Sharp 


20.0 





45.0 


5.0 


13.5 


70 


90 


96 


5.0 


.5 



(D 
(A 

n 



c/) 
n 



(A 

<D 



O 
(A 



Table 2. Summary of pier-scour measurements — Continued 



BSDMS^ 
site 


Measure 
ment 
no.2 


Date 


Pier 
identi- 
fication 


Scour 
con- 
dition 


Pier 
shape 


Pier 

L 
(ft) 


Pier skew 
(deg) 




Pier 
b 
(ft) 


V 

(ft/s) 


Flow depth 
(ft) - 


Particle size for which 

indicated percentage of 

bed material isfiner^ 


Scour 

depth 

(ft) 


Meas- 
ure- 
ment 
error 
(ft) 




3 


-0 
CO 




no. 


D50 
(mm) 


D90 

(mm) 


D95 

(mm) 
















Alas 


;ka 


-continued 

















c 


1 


8 


08/11/71 


P4 


Clear 


Sharp 


20.0 









"5.0 


9.5 


17.5 


70 


90 


96 


5.0 


1.0 


m 


4 


29 


09/02/71 


PI 


Live 


Round 


- 









^15 


9.5 


12.0 


90 


130 


140 


5.0 


.5 




4 


30 


09/04/71 


PI 


Live 


Round 


- 









^15 


11.5 


15.0 


90 


130 


140 


5.5 


.5 





















IVIaryli 


and 
















3" 


23 


153 


07/13/90 


Left 


Clear* 


Sharp 


41.7 









5.0 


7.7 


7.9 


110 


290 


350 


1.1 


1.0 







23 


154 


04/01/93 


Left 


Clear* 


Sharp 


41.7 









5.0 


6.8 


6.8 


110 


290 


350 


1.4 


1.0 


ai 

CA 


23 


155 


07/13/90 


Right 


Clear* 


Sharp 


41.7 









5.0 


8.6 


9.9 


110 


290 


350 


2.7 


1.0 


? 


23 


156 


04/01/93 


Right 


Clear* 


Sharp 


41.7 







Ohi( 


5.0 
3 


6.2 


8.0 


110 


290 


350 


1.7 


1.0 


00 
<D 

a. 


44 


288 


12/19/90 


P2 


Clear* 


Round 


24.3 









2.5 


3.7 


5.6 


60 


74 


75 


.7 


.3 


CD 

ai 

3 

(A 


















Virginia 
















53 


358 


05/03/89 


P2 


Clear* 


Round 


41.0 









2.0 


1.8 


1.5 


72 


220 


250 


.8 


1.0 




53 


359 


05/07/89 


P2 


Live* 


Round 


41.0 









2.0 


5.1 


2.2 


72 


220 


250 


.6 


1.0 




53 


360 


04/22/92 


P2 


Live* 


Round 


41.0 









2.0 


5.2 


5.5 


72 


220 


250 


1.6 


1.0 




53 


361 


05/03/89 


P3 


Live* 


Round 


41.0 









2.0 


4.0 


4.0 


72 


220 


250 


1.0 


1.0 




53 


362 


05/07/89 


P3 


Live* 


Round 


41.0 









2.0 


5.3 


5.0 


72 


220 


250 


1.2 


1.0 




53 


363 


04/22/92 


P3 


Live* 


Round 


41.0 









2.0 


8.5 


8.6 


72 


220 


250 


2.5 


1.0 




55 


376 


03/29/91 


P2 


Clear 


Round 


30.0 









2.0 


3.7 


2.5 


55 


95 


110 


1.5 


1.0 




55 


377 


06/05/92 


P2 


Clear 


Round 


30.0 









2.0 


5.5 


10.5 


55 


95 


110 


2.1 


1.0 




55 


378 


03/24/93 


P2 


Clear 


Round 


30.0 









2.0 


6.4 


10.5 


55 


95 


110 


1.8 


1.0 





Data compiled from U.S. Geological Survey, Bridge Scour Data Management System (BSDMS) or from Landers and Mueller (1996, p. 58-66). 

Measurement number in Landers and Mueller (1996). 

Where Dgg or Dg^ were not available, they were calcul ^^ 

Substantial debris on pier could lead to underestimation of scour depths. 

Maximum pier width in Landers and Mueller (1996). 
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HEC-18/Mueller equation development. Therefore, a substan- 
tial portion of the data used in this study could be considered 
independent for evaluating the HEC-18/Mueller equation. The 
overlap between data used in this evaluation and in development 
of each of the other four equations is unknown, but is probably 
smaller because the other equations are less recent. In addition, 
this study focused on the subset of pier-scour data pertaining to 
coarse-bed streams. Data used for evaluation of the equations 
included only data from streams where the bed-material D^q 
was greater than 50 mm, with an average D^q equal to 76 mm. 
In contrast, the data used to develop the HEC-18/Mueller equa- 
tion included Dgg values as small as 0.18 mm, with an average 
D50 equal to 12 mm. 

Few pier-scour data are available for rare floods in coarse- 
bed streams. Data used in this study were collected during 
floods substantially smaller than the 100-year event, with the 
notableexceptionof the Yellowstone River (table 1). However, 
predicted scour depths for the 100-year and even 500-year 
floods typically are required for bridge design. More scour 
measurements from flows closer to the 100-year and 500-year 
events are needed in order to test the performance of the scour- 
prediction equations for these less common events. 



Description of Pier-Scour Equations for 
Coarse-Bed Streams 



The five pier-scour equations evaluated in this report 
include the Simplified Chinese equation, the Froehlich equa- 
tion, the Froehlich design equation, the HEC-18/Jones equa- 
tion, and the HEC-18/Mueller equation. All of these equations 
attempt to account for the effects of the bed-material size on 
scour. The equations have been modified from their original 
format so that all units are in the foot-pound-second or English 
unit convention. 



b is the width of bridge pier, in feet; 

yg is the depth of flow directly upstream from the pier, in 
feet; 

D„ is the mean particle size of the bed material, in feet 
(for this study D^q was used for D^); 

Vg is the approach velocity directly upstream from the 
pier, in feet per second; 

Vg is the critical (incipient motion) velocity, in feet per 
second, for the D„-sized particle. If the density of 
water is assumed to be 62.4 pounds per cubic foot, and 
the bed material is assumed to have a specific gravity 
of 2.65, the equation for Vg can be expressed as: 



V = 3.28 



8.85£)„ +6.05E 



■ 10 H- 0.3048);„ ■ 
(0.3048A„ )°-^2 



0.5 



(la) 



Vjg is the approach velocity, in feet per second, corre- 
sponding to critical velocity at the pier. V,-^ can be cal- 
culated using the following equation: 



/£) \ 0-053 

y,., = 0.645 [-^J K 



(lb) 



The Simplified Chinese Equation for live-bed scour (when 
Vg>Vg) is expressed as: 



r, ocn I' 1,0.6 0.15 -0.07 fYo-}jc'\'' 
y, = 0.950K,b y^ D,„ \ _ 

c ic 



(2) 



where the exponent c is calculated using the following equation: 



y. 8.20 + 2.23 logD^ 



V, 



(2a) 



Simplified Chinese Equation 

The Simplified Chinese pier-scour equation is based on 
laboratory and field data from China (Landers and Mueller, 
1996, p. 98-100). This equation has different forms depending 
upon whether the scour condition is live-bed scour (bed mate- 
rial upstream from bridge is in motion) or clear-water scour 
(bed material upstream from bridge is not in motion). The Sim- 
plified Chinese equation for clear-water pier scour is defined as: 



J. 



11/111' J,0.6 0.15 -0.07/' V, 

1. 141 K^^b y„ D^ 



K 



y -y. 



(1) 



where: 



is the depth of pier scour below ambient bed, in feet; 
is the simplified pier-shape coefficient: 
Kg = 1.0 for cylinders, 

= 0.8 for round-nosed piers, 

= 0.66 for sharp-nosed piers; 



To determine whether the equation for clear-water scour 
(equation 1) or live-bed scour (equation 2) was appropriate, crit- 
ical velocity (Vg) was calculated using equation la. Pier scour 
was calculated using equation 1 for measurements where Vg<Vg 
and equation 2 for measurements where Vg>Vg, regardless of 
the scour condition reported in the BSDMS. 

Froehlich Equation 

The Froehlich equation (Froehlich, 1988) was derived 
using regression analyses of pier-scour data from several inves- 
tigations (Landers and Mueller, 1996, p. 101) and is defined as: 



/ i,n0.62 /,, nO.46 



0.08 



(3) 



where: 

yg, yg, and b are as previously defined; 
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<!> is a dimensionless coefficient based on the shape of 
the pier nose, as follows: 
(/) = 1.3 for square nosed-piers, 
= 1.0 for round-nosed piers, 
= 0.7 for sharp-nosed piers; 
Fri is the Froude Number directly upstream from the pier, 

dimensionless; 
bg is the width of the bridge pier projected normal to the 

approach flow, in feet; and 
D50 is the particle size for which 50 percent of the bed 
material is finer, in feet. 
According to Landers and Mueller (1996, p. 101) the 
Froehlich equation is based on field measurements that were 
presumed to have been made under live-bed scour conditions. 
However, in this study, the equation also was applied at all of 
the sites, even though clear-water scour was reported for many 
measurements. 



K2 is the correction factor for angle of attack of flow 

at the pier, dimensionless; 

Kj, is the correction factor for bed condition, dimen- 

sionless; 

^4 is the correction factor for armoring of coarse- 

bed material, dimensionless; and 

b and Fri are as previously defined. 

The Ki, K2, and Kt, correction factors are defined in table 3. ^4 
is calculated by the following equation: 



K,= [1 -0.89(1 -y«)']°' 



where: 

Vjf is the velocity ratio, dimensionless, calculated as: 



(6) 



Froehlich Design Equation 



v„ 



r y. - y,'= 



c90 



V, 



;50- 



(6a) 



A scour-depth estimation method based on regression 
analysis, where underestimates are as likely as overestimates, is 
undesirable for bridge design because underestimation of scour 
depth is not acceptable. Froehlich (1988) found that by adding 
pier width (b) to the scour depth computed by equation 3, pier 
scour was not underestimated for any of the bridges in the data 
set used. Thus, the Froehlich design equation (Mueller, 1996, 
p. 102) is defined as: 



Vg is as previously defined; 

y,-50 is the approach velocity, in feet per second, required 
to initiate scour at the pier for the particle size D^q. 
Vi^Q is calculated as follows: 



y,5o=0.645r^| V. 



c50 



(6b) 



3;, = 032b^Fr°-H^ 



0.62 



0.46 



a 



50 



+ b 



(4) 



where all variables are as previously defined. The Froehlich 
design equation is included as a pier-scour calculation option 
within the computer model HEC-RAS, Version 3.1 (Brunner, 

2002). 



D50 is as previously defined, 

y^so is the critical velocity, in feet per second, for incipi- 
ent motion of the particle size D50, and is further 
defined as follows: 



y,3o= ii.2iy>5o^' 



(6c) 



HEC-1Q/Jones Equation 

The HEC-18/Jones equation is based on the Colorado 
State University (CSU) equation (Richardson and others, 1993). 
When compared with data from field measurements (Richard- 
son and others, 1993), the CSU equation was found to more reli- 
ably calculate pier scour when compared to several other equa- 
tions. The HEC-18/Jones equation incorporates a correction 
factor K^ to account for armoring of the scour hole. Richardson 
and Davis (1995, p. 36-38) define the HEC-18/Jones equation 



y, = 2.0y „K,K,K,K,[^ 



0.65 



Fr 



0.43 



(5) 



where: 

yg and y^ are as previously defined; 
Ki is the correction factor for pier-nose shape, 

dimensionless; 



y^9Q is the critical velocity, in feet per second, for incipi- 
ent motion of the particle size Dgg, and is 
calculated as follows: 



y,9o= 11.21y„'^X''' 



(6d) 



D90 is the particle size for which 90 percent of the bed 
material is finer, in feet. 

Mueller and Wagner (in press) report that at sites where 
only D50 and D^ are determined, DgQ can be calculated by: 



D -D ("^ 

^90 - ^50^^^ 



n \ 1-282 



(6e) 



where the ratio — ^ is referred to as the gradation coefficient 
(Mueller, 1996). '° 
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Table 3. Ki, /C2, and K^ correction factors for the HEC-18/Jones equation 

[Information in table from Richardson and Davis (1995). Abbreviations: b, pier width, in feet; 
ft, feet; L, pier length, in feet. Symbols: >, greater than; >, greater than or equal to] 



Correction Factor 



Value or Equation 



K^: Pier-nose shape 

Square nose 
Round nose 
Circular cylinder 
Group of cylinders 
Sharp nose 



1.1 
1.0 
1.0 
1.0 
.9 



K2: Angle of attack of flow (G) 

^3: Bed condition 

Clear-water scour 

Live-bed scour: 

Plane-bed and antidune bedform 
Small dunes, 3 ft > dune height > 0.6 ft 
Medium dunes, 9 ft > dune height > 3 ft 
Large dunes, dune height > 9 f t 



[CosO + (L/b) Sine]""\ 

where 6 = angle of attack of flow 

If L/fo>12, use 12 as a maximum 



l.I 

l.I 

I.I 

l.I to 1.2 

1.3 



Equation 6 is based on research indicating that when the 
approach velocity (Vg) is too low to move the D90 size of the 
bed material, scour depth is reduced (Richardson and Davis, 
1995, p. 37). In this situation {¥„< V^gQ, equation 6a), K4 will 
be less than 1.0. However, for comparison purposes in this 
study, 7^4 was allowed to be less than 1.0 even when Vg> Vj^q. 
Richardson and Davis (1995, p. 38) recommend a minimum 
value of 0.7 for the HEC18/Jones A'4 correction factor and indi- 
cate that when the velocity ratio (V^) exceeds 1.0, K/^ should 
default to 1.0. Richardson and Davis (1995, p. 38) also suggest 
that the HEC-18/Jones equation be applied only to sites where 
D50 is greater than 60 mm. However, the equation also was 
used in this study for four scour measurements from sites 55 and 
85 where the bed-material Dgg was estimated to be between 50 
mm and 60 mm. 



the D^Q sized particles (V^so) that is slightly different from the 
y^50 used in the HEC-18/Jones equation. Mueller's equation 
for y^jQ varies in format for different Dgg sizes (Mueller, 1996, 
p. 160). For Dgo greater than 20 mm, V^^q is calculated as fol- 
lows: 



y,3o = 12.21yr05o' 



(7a) 



where all terms are as previously defined. 

The y,-95 term is calculated like Vf^Q in equation 6b as fol- 
lows: 



.£) s 0.053 

y,„ = 0.645 (-^J V, 



(7b) 



HEC-18/Mueller Equation 

Mueller (1996, p. 160) proposed a modified K^ correction 
factor for the HEC-18 equation: 



K^ = 0.4 



V„-V, 



!50 



V,. 



V, 



0.15 



(7) 



where V/gj is the approach velocity, in feet per second, required 
to initiate scour at the pier for the Dg^ particle size, in feet, and 
all other terms are as previously defined. Mueller (1996, p. 160) 
suggests an equation for critical velocity for incipient motion of 



V^gs, is calculated similarly to V^gg in equation 7a as follows: 

y,„5 = 12.21V ^'<' (7c) 



where Dg^ is the particle size for which 95 percent of the bed 
material is finer, in feet. 

Mueller and Wagner (in press) determined that when 
information about Dg^ is not available, Dg^ can be calculated 
by: 



Do 






84 



1.645 



(7d) 
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Richardson and Davis (2001, p. 6.6) recommend a mini- 
mum K/^ value of 0.4. Furthermore, J. Sterling Jones (Federal 
Highway Administration, oral commun., 2002) recommends 
that measured scour be compared to calculated scour only if the 
measured approach velocity (Vg) exceeds 0.4y;5o because bed 
material theoretically would not be scoured at smaller approach 
velocities. Scour measurements made for approach velocities 
less than 0.4y;5Q likely reflect scour holes that are remnants of 
an earlier flood. For this study, one measurement where 
Vg<0.4Vj^Q is included for comparison purposes (site 53, mea- 
surement 358, table 2). 

When the approach velocity exceeds the critical velocity 
for the coarsest size fraction (Vo>V,c95 or V,c9o)' J- Sterling 
Jones (Federal Highway Administration, oral commun., 2002) 
recommends that K^ default to 1.0. However, the HEC-18/ 
Mueller equation for K/^ (equation 7) is based on data that 
include scour measurements where Vo> ^,^95 (David S. Mueller, 
U.S. Geological Survey, oral commun., 2002). Therefore, 
when comparing measured scour to calculated scour in this 
study, K/^ was calculated by equation 7, even when Vg>Vjggy 
However, when designing or analyzing a bridge, K/^ values 
should typically be no lower than 1.0 when Vg exceeds V,(,95- 

Table 4 summarizes the characteristics of the data used to 
develop equation 7 (Mueller, 1996, p. 93). The HEC-18/Muel- 
ler equation might not be reliable when applied outside the 
ranges of variables listed in table 4. Furthermore, according to 
the Federal Highway Administration (J. Sterling Jones, oral 
commun., 2002) certain Djg/Dgg combinations can provide 
erratic results when using equation 7, and further refinement of 
the K/^ correction factor is being evaluated. The HEC-18/Muel- 
ler equation is included as a pier-scour calculation option within 
the computer model HEC-RAS, Version 3.1 (Brunner, 2002). 



Evaluation of Pier-Scour Equations for 
Coarse-Bed Streams 



sured values, but plot close to the line labeled "calculated pier 
scour equals measured pier scour." Similarly, the boxplots of 
residuals (calculated scour minus measured scour) in figure 3 
can be used to visually assess reliability of the five equations. 
On this basis, the most reliable equations are associated with 
boxes that are short, but that generally lie above the line indicat- 
ing zero residual depth of scour. 

Some error typically is involved in the measurement of 
pier scour; this error usually is estimated at the time the mea- 
surement is made. For the coarse-bed streams in this study, 
measurement errors ranged from 0.3 ft to 1.0 ft. In some 
instances, the measurement error was greater than the calcu- 
lated residual (tables 2 and 5). 



Simplified Chinese Equation 

Pier-scour depths calculated using the Simplified Chinese 
equation were smaller than measured scour depths for 15 of the 
18 measurements from Montana and for 14 of the 24 measure- 
ments from the 4 other States (table 5, figs. 2A and 3). How- 
ever, 86 percent of the underestimated scour depths came within 
2 ft of the measured scour. This equation occasionally pro- 
duced pier-scour estimates that were less than zero. Calculated 
negative pier-scour values are reported in table 5 and plotted in 
figure 2 A as zero. Calculated negative pier scour values also 
are assumed to be zero for purposes of calculating residuals 
reported in table 5 and shown in figure 3. Calculated scour 
depth was less than measured scour depth by about 2 to 5 ft at 
P4, site 1 in Alaska, where substantial debris was reported for 
measurements 7 and 8. Debris on the pier probably contributed 
to the measured scour and accounted for the underestimation by 
the equation. For the 13 measurements where calculated scour 
was greater than measured scour, the difference between calcu- 
lated and measured scour ranged from 0.2 to 1.6 ft. 



Pier-scour depths calculated by the five equations were 
compared to 42 pier-scour measurements. In this section, com- 
parisons for the results from each equation are discussed and 
presented using scatter and boxplots. 

The equations were evaluated using two criteria, both of 
which are important to bridge design and to analysis of existing 
structures: 

1 . Number and magnitude of underestimates need to be min- 
imal, and 

2. Calculated scour depths need to match measured scour 
depths as closely as possible. 

Results indicating the performance of the five equations 
are presented in table 5, and illustrated in figures 2 and 3. The 
scatterplots in figure 2 compare calculated pier scour for each 
of the five equations to measured pier scour for each measure- 
ment at each site. Generally, the most reliable equations result 
in calculated pier scour values that are mostly greater than mea- 



Froehlich Equation 

Results from the Froehlich equation were similar to those 
from the Simplified Chinese equation (table 5, figs. 2A, 2B, and 
3). The Froehlich equation is based on a regression analysis, 
where underestimates are as likely as overestimates. Pier-scour 
depths calculated using the Froehlich equation were smaller 
than measured scour depths for 14 of the 18 measurements from 
Montana and for 15 of the 24 measurements from the other 
States. Seventy-six percent of the underestimates were within 
2 ft of measured scour. 

Substantial debris on the pier at site 1 in Alaska likely 
explains why calculated scour depth was substantially less than 
measured scour depth for measurements 7 and 8. For the 12 
measurements where calculated scour was greater than mea- 
sured scour, the difference between measured and calculated 
scour ranged from 0.1 to 0.6 ft. 
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Table 4. Summary of scour measurement data used to develop the HEC-18/Mueller equation 

[Dj is the particle size for which "x" percent of the bed material is finer. V^jg '* the critical incipient motion velocity for particles of size 
D50; >'o depth of flow directly upstream from the pier, in feet. Abbreviations: ft, foot; ft/s, foot per second; mm, millimeter; mi , square mile] 



Variable 


Units 


Number of 
measure- 
ments 


Minimum 


Median 


Maximum 


Mean 


Standard 
deviation 


Depth of scour, Vj 


ft 


224 





2.0 


25.1 


3.1 


3.5 


Approach velocity, Vg 


ft/s 


224 


.6 


3.7 


14.7 


4.6 


3.1 


Approach depth, >•„ 


ft 


224 





13.4 


39.2 


14.5 


8.9 


O50 


mm 


224 


.18 


.84 


108 


12 


23 


Ds4 


mm 


220 


.37 


3.3 


233 


27 


43 


Dgs 


mm 


195 


.48 


7.0 


350 


33 


65 


Gradation coefficient 


dimensionless 


220 


1.3 


2.3 


12 


3.5 


2.5 


Pier width, b 


ft 


224 


1.2 


3.4 


15.1 


4.4 


2.8 


Drainage area 


mi2 


156 


64 


1,420 


60,700 


7,670 


16,500 


Channel slope 


dimensionless 


178 


.00016 


.00075 


.0050 


.0010 


.0011 


Jo/b 


dimensionless 


224 


.12 


3.4 


14 


4.0 


2.9 


WO50 


dimensionless 


224 


8.5 


990 


14,200 


1,890 


3,180 


Froude number 


dimensionless 


224 


.039 


.22 


.83 


.24 


.17 


Vo/Vc50 


dimensionless 


224 


.44 


1.1 


5.4 


1.3 


.9 


yjb 


dimensionless 


224 





.6 


2.1 


.7 


.4 



'Data from Mueller (1996). 
"Gradation coefficient = D^JD^, 



The Froehlich equation calculated about 1 ft of scour for 
all the floods at the Montana sites, even though the measured 
scour varied from 0.3 to 5.5 ft. The calculated scour most likely 
did not vary because the Montana bridges have similar bed- 
material sizes and pier widths (table 2). 



tion. The HEC-18/Jones equation overestimated pier scour at 
site 4 in Alaska by 14.9 and 18.8 ft (measurements 29 and 30, 
respectively, table 5). As discussed above, this large overesti- 
mation is likely due to the 15-ft-wide pier at this bridge. Over- 
estimates at other sites ranged from 0.3 to 9.2 ft. 



Froehlich Design Equation 

Pier-scour depths calculated using the Froehlich design 
equation exceeded measured scour depths for every measure- 
ment except measurements 231 and 232, at site 32 in Montana 
(table 5, figs. 2C, and 3). Pier scour calculated using the 
Froehlich design equation for measurement 23 1 was only 0.6 ft 
less than measured pier scour, and for measurement 232, calcu- 
lated scour was equal to measured scour. Pier-scour estimates 
obtained by this equation tend to be greatly influenced by the 
pier-width variable (fig. 2C). 

At site 4 in Alaska, pier-scour depths calculated using the 
Froehlich design equation were more than 15 ft greater than 
measured scour depths. This large overestimation is likely due 
to the 15-ft-wide pier. Overestimates from the Froehlich design 
equation at all other sites ranged from 0.1 to 5.4 ft. Overesti- 
mates calculated using the HEC-18/Jones equation and the 
HEC-18/Mueller equation also were larger at site 4 than at all 
other sites. 



HEC-18/Jones Equation 

Pier-scour depths calculated using the HEC-18/Jones 
equation exceeded measured pier scour for all 42 observations 
(table 5, figs. 2D and 3). Many overestimates were larger for 
the HEC-18/Jones equation than for the Froehlich design equa- 



HEC-1Q/Mueller Equation 

Pier-scour depths calculated using the HEC-18/Mueller 
equation were seldom less than measured pier-scour depths 
(table 5, figs. 2E and 3). Calculated scour depths were less than 
measured scour depths for 8 of the 42 observations, but 2 of 
these 8 underestimates were probably due to substantial debris 
on the pier at site 1 (measurements 7 and 8) in Alaska. For the 
other six observations where scour was underestimated, the dif- 
ference between calculated and measured scour ranged from 0.1 
to 1.1ft. 

The HEC-18/Mueller equation resulted in lower K^ values 
than did the HEC-18/Jones equation, and the minimum allow- 
able K^ value (0.4) in the Mueller version is lower than the min- 
imum K^ value (0.7) recommended by the FHWA for the HEC- 
18/Jones equation. Therefore, the HEC-18/Mueller equation 
resulted in smaller calculated pier-scour depths than the HEC- 
18/Jones equation. Calculated scour depths from the HEC-18/ 
Mueller equation were in better agreement with measured scour 
depths compared to calculated scour depths by the other equa- 
tions that did not underestimate pier scour. 

At site 4 in Alaska, where the pier width is 15 ft, the HEC- 
18/Mueller equation overestimated pier scour by 3.9 and 4.5 ft 
(tables 2 and 5). This disparity might be explained by the fact 
that 15 ft is the maximum pier width in the data set used to 
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Figure 2. Comparison of calculated to measured pier scour for five equations. Data from U.S. Geological Survey 
Bridge Scour Data Management System (Chad R. Wagner, U.S. Geological Survey, written commun., 2002). 
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Figure 3. Distribution of residual depth of pier scour (calculated scour minus the measured scour). Residuals 
less than 0.0 indicate pier scour was underestimated. 



Table 5. Comparison of calculated pier scour to measured pier scour from five equations 

[All data for measured and calculated pier scour, measurement error, and residual are reported in feet. Pier identification: number or location of pier at which measurement took place. Residual: 
calculated scour minus measured scour, negative residual values indicate underestimated scour. Figures might not add to totals because of independent rounding. Abbreviations: BSDMS, Bridge 
Scour Data Management System. Symbol: ~, not applicable] 
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Table 5. Comparison of calculated pier scourto measured pier scourfrom five equations — Continued 
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Measurement number from Landers and Mueller (1996). 

Simplified Chinese equation resulted in negative values for calculated pier scour; calculated scour is reported as zero. 
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develop the HEC-18/Mueller equation (table 4). As discussed 
above, the Froehlich design equation and the HEC-18/Jones 
equation overestimated pier scour at site 4 by 14.9 to 18.8 ft. 
Richardson and Davis (2001, p. 6.7) define a correction factor, 
K^^;, that adjusts the HEC-18/Mueller equation for wide piers. 
However, the ratio of pier width to the Dgg-particle size (b/D^o) 
at site 4 did not meet the criteria for using the K^^, correction fac- 
tor. In addition, Richardson and Davis recommend that K^^, be 
used with caution because K^^, is based only on few data from 
flume experiments. For these reasons, K^^, was not applied to the 
calculations in this study. 



The HEC-18/Mueller equation resulted in relatively few 
underestimates, with an average underestimate of 0.9 ft. For 
measurements where the HEC-18/Mueller equation overesti- 
mated scour, the average overestimate was 1 .4 ft. This average 
overestimate was smaller than overestimates from the other 
equations that did not underestimate scour (Froehlich design 
and HEC-18/Jones equations). 



Summary and Conclusions 



Summary Statistics for Calculated and 
Measured Pier Scour 



Statistics for calculated and measured pier scour are sum- 
marized in table 6. The averages of pier-scour depths calculated 
from all five equations ranged from 1.5 to 6.1 ft; the average of 
measured pier-scour depths was 2.2 ft. Average calculated pier- 
scour depths from the Froehlich equation (1.5 ft). Simplified 
Chinese equation (1.7 ft), and the HEC-18/Mueller equation 
(3.1 ft) were relatively close to the average of measured pier- 
scour depths (2.2 ft). Although averages of calculated pier- 
scour depths from the Froehlich and Simplified Chinese equa- 
tions were closer to the average of measured pier scour depths, 
these equations resulted in the most underestimates of all the 
equations. 



Five pier-scour equations were evaluated based on existing 
pier-scour, bed-material, bridge, and streamflow data for 
coarse-bed streams in Montana, Alaska, Maryland, Ohio, and 
Virginia. Pier scour-depths calculated for flood events with 
approximate recurrence intervals of less than 2 to 100 years by 
the Simplified Chinese equation, the Froehlich equation, the 
Froehlich design equation, the HEC-18/Jones equation, and the 
HEC-18/Mueller equation were compared to 18 pier-scour 
measurements at 3 bridge sites in Montana and 24 pier-scour 
measurements at 6 bridge sites in 4 other states. Site informa- 
tion and measurements used in this study are summarized in the 
report. 

When applied to data from the bridge sites, the equations 
produced results that were consistent between Montana and the 
four other states. The Simplified Chinese and Froehlich 



Table 6. Summary statistics for calculated pier scour for five pier-scour equations and measured pier scour 

[Statistics include data for all 42 sites. Symbol: ~, not applicable] 
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1.2 
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Froehlich design'* 


5.6 
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.6 
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3.5 
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HEC-18/Jones 


6.1 


4.4 
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3.9 


42 


HEC-18/Mueller^ 


3.1 


1.8 


.9 


8 




1.4 


32 


Measured 


2.2 


1.5 


- 


- 




- 
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Data from table 5. 

Absolute value of average of differences between scour calculated using indicated equation and measured scour for measurements where scour was underestimated. 
Average of differences between scour calculated using indicated equation and measured scour for measurements where scour was overestimated. 
^he Simplified Chinese, Froehlich, and Froehlich design equations, each accurately estimated scour for one measurement (residual values were equal to 0.0, table 5). 
The HEC-1 8/Mueller equation accurately estimated pier scour for two measurements (residual values were equal to 0.0, table 5). 
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equations both underestimated scour at several sites. Eighty-six 
percent of the underestimates from the Simplified Chinese 
equation were less than 2 ft, and 76 percent of the underesti- 
mates from the Froehlich equation were less than 2 ft. The 
HEC- 1 8/Jones equation resulted in the largest overestimates for 
the equations examined in the study. The Froehlich design 
equation also overestimated scour, although overestimates were 
smaller than for the HEC-18/Jones equation. 

The HEC-18/Mueller equation generally predicted scour 
better than the other four equations according to the evaluation 
criteria used in this study. Scour was seldom underestimated, 
and calculated scour was closer to measured scour than for the 
other equations that did not underestimate scour. This equation 
might not be reliable when applied to sites where variables are 
outside the ranges within the dataset from which the equation 
was derived. This equation is being revised and sometimes can 
produce erratic results for certain combinations of D^Q/Dg^. 
Measurements of pier scour at additional sites are needed to fur- 
ther evaluate the HEC-18/Mueller equation. Additional scour 
measurements during rarer events, such as the 100-year and 
500-year flood, also are needed to further test the equation. 
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